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Genetic learning algorithms are widely used to control ultrafast optical pulse shapes for photo- 
induced quantum control of atoms and molecules. An unresolved issue is how to use the solutions 
found by these algorithms to learn about the system's quantum dynamics. We propose a simple 
method based on covariance analysis of the control space, which can reveal the degrees of freedom in 
the effective control Hamiltonian. We have applied this technique to stimulated Raman scattering 
in liquid methanol. A simple model of two-mode stimulated Raman scattering is consistent with 
the results. 
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The central challenge of coherent control of quantum 
dynamics is to find the optimal path to guide a quantum 
system from its initial state to some target final state |l|, 
0. Several theoretical methods have been developed to 
aid this search 0, 0, and there has been considerable 
experimental success as well|5j,|y|. However, in all but the 
simplest systems, the search is hampered by incomplete 
knowledge of the system Hamiltonian. Strongly coupled 
systems such as large molecules in condensed phase are 
so complicated that it is nearly impossible to calculate 
optimal pulse shapes in advance. 

Feedback learning algorithms overcome this limitation 
by using the physical system itself to explore its own 
quantum dynamics through an experimental search 0. 
A typical search experiment compares the ability of sev- 
eral thousand different shaped laser pulses to transform 
the system from its initial state at time t = to some 
desired target state \\) at a later target time t ~T. Ex- 
amples of transformations that have been studied include 
molecular photodissociation, atomic photoexcitation and 
photoionization. The pulse shapes are selected through 
a fitness-directed search protocol, such as a genetic al- 
gorithm The fitness is a measured quantity propor- 
tional to the objective functional J[H;xi], which is the 
square of the projection of x) onto at the end of the 
experiment: 

J[H;x i ] = \( X (T)\i,(T))\ 2 . (1) 

J depends on the Hamiltonian H for the system evolu- 
tion, which depends in turn on the laser electric field E(t) 
determined by the settings of the n pulse shape control 
parameters Xi, i = 1 . . .n. J reaches its extreme value 
for the optimal pulse. This pulse can be calculated using 
optimal control theory if H is known 0; otherwise, it 
must be discovered through the learning search. 

Several recent papers have suggested modifications or 
extensions of learning feedback that can measure prop- 
erties of the system Hamiltonian 0, 0, [nl [^. Here 
we propose a different approach based on analysis of the 
trial experiments. We will show that the ensemble of trial 
pulse shapes can reveal essential features of the dynam- 
ics. 



Genetic algorithms and similar evolutionary search 
strategies have many different variations [l3l | . Our imple- 
mentation starts with approximately 50 randomly gener- 
ated optical pulse shapes produced by spectrally filter- 
ing an ultrafast laser pulse l4]. Each pulse shape is 
described by a column matrix of control parameters Xi 
called a genome consisting of about 25 numbers {genes), 
each encoding the amplitude and/or phase of a differ- 
ent segment of the optical spectrum. The control target 
is measured for each pulse shape. Then the algorithm 
creates a new generation of pulse shapes by combining 
attributes of the fittest members of the previous gener- 
ation 15]. After several generations, the pulse shapes 
usually cluster near high fitness regions of the search 
space. When the algorithm finds a pulse shape or several 
shapes that cannot be improved over many generations, 
the search stops, and the highest fitness pulse shape is de- 
clared the optimal solution to the search. We test 1000 to 
10,000 pulse shapes in a typical experiment. We main- 
tain a record of every pulse shape, its fitness, and its 
parentage (genealogy). 

The learning algorithm achieves control without prior 
knowledge of the system Hamiltonian, and has far more 
degrees of freedom n than the minimum required for con- 
trol. The number of possible solutions is exponential in 
n. In a typical search, the phase of each color is adjusted 
by the spectral phase filter to a precision of about ten de- 
grees, so there are 2 5 possible values of each gene. This 
means that the number of possible solutions for a genome 
of length 25 is 2 5x25 ~ 4 x 10 37 . Genetic algorithms 
can search this large state space with great efficiency |8|. 
Unfortunately, simply finding a good solution has not of- 
ten provided significant insight into the system dynamics 
or Hamiltonian. The optimal pulse shape found by the 
learning algorithm, while sufficient to achieve control, is 
often complicated and may contain unnecessary features. 

The conditions for reaching an extremum in J[H; Xi] 
may only depend on two or three essential features of the 
control field E(t). These features are not obvious in the 
successful genome because they may depend on all 25 
genes. The Hamiltonian could be written in a much sim- 
pler form if these essential degrees of freedom (v,j ) were 
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found. Here we show how to establish the Uj through 
covariance analysis of the pulse shapes evaluated during 
the learning search^^ . Covariance analysis is commonly 
used to reduce the dimensionality of and to find patterns 
in high dimensional data sets. 

We propose to apply these techniques, not to sets of 
data, but to the control space for the experiment. Lin- 
ear combinations of genes with high fitness should ap- 
pear correlated in the fitness-driven genetic algorithm. 
These correlated linear combinations correspond to the 
principal components of the control field that direct the 
quantum dynamics under investigation. Principal con- 
trol analysis is the application of covariance techniques 
to a fitness directed search. 

Principal control analysis is implemented on our sys- 
tem by calculating the covariance matrix of the set of all 
pulse shapes in the search, defined by: 



sential features of the search solutions. The objective 
functional also takes on a simpler form in this basis: 



C ij = (6 i 5 j )-{8 i )(S j ), 



(2) 



where the expressions Si = a?i+i —Xi,i = l...n—1 are the 
nearest neighbor phase differences. By using the phase 
differences in the analysis we remove the ambiguity as- 
sociated with the unimportant global phase. The covari- 
ance matrix is not the only measure of correlation. We 
could also weight the terms of the covariance matrix us- 
ing the fitness, or normalize each term to the individual 
gene variance. In this paper, we will use the simple co- 
variance. This is appropriate because all the S^s are of 
the same type. 

Once the covariance matrix is determined, we calcu- 
late its eigenvectors and eigenvalues. Each eigenvalue Ad- 
measures the variance of the projections onto the cor- 
responding eigenvector. This has a special meaning for 
a learning control search: it shows how far the control 
setting moved during the learning process. A small sub- 
set of eigenvalues usually contains most of the weight of 
the trace of the covariance matrix. The correlations of 
the projections with the pulse shape fitnesses allow us 
to determine which control directions were most impor- 
tant for the physical process under consideration. The 
controls expressed in the basis of the eigenvectors are 
uncorrelated: Each of these controls changes the fitness 
without correlation with the others over the search set. 

We propose that these eigenvectors with the largest fit- 
ness correlation are the essential control directions (uj). 
We expect the eigenvectors with the larger eigenvalues 
to be most strongly correlated with the fitness of a pulse 
shape solution. Conversely, a low correlation indicates 
those eigenvectors that have not contributed substan- 
tially to increasing the fitness during the search. These 
eigenvectors correspond to extraneous dimensions, which 
could be eliminated (i.e., their projection set to zero) 
without losing substantial control. 

By projecting the GA solutions onto the k < n eigen- 
vectors that correlate best with the fitness (the principal 
controls), we reduce the dimension of the control space. 
The solutions with highest fitness, when expressed in the 
reduced basis of the principal controls, represent the es- 



J = J[H;ui t ... t k] 



(3) 



where each Uj can assume a range of values on the order 
of ±WAj, The search for the specific target state is now 
a matter of optimizing each control Uj over this range. 

In summary, we propose to apply covariance tech- 
niques to the control space of learning feedback exper- 
iments. Control values derived from the genomes are 
analyzed by a covariance matrix, defined in Eq. |5J The 
matrix eigenvectors are independent control directions. 
The correlation of the fitness with the eigenvectors sug- 
gests which control directions are the most important. 
The corresponding eigenvalues indicate the necessary ex- 
cursion along each axis. Therefore, searches conducted in 
the eigenvector basis are more efficient. Finally, the best 
solutions are found by projecting the optimal learning 
control solutions onto the important control directions. 
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FIG. 1: Left: Raman spectra (a) and Wigner representation 
(c) of the optimal pulse shape after optimizing the symmetric 
C-H stretch mode. Wigner representation (e) of the control 
pulse shape found by principal control analysis. Right: Simi- 
lar for antisymmetric stretch. 

We now apply this analysis to a well-studied control 
problem: the selective excitation of vibrational modes 
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in liquid methanol. The experiment has been described 
previously 0, 0, 0|. An intense shaped 800 nm ul- 
trafast laser pulse (the pump laser) is focused into a 
cell containing methanol. Above a threshold fluence, 
the pump induces stimulated Raman scattering into ei- 
ther the symmetric or antisymmetric Raman-active C-H 
stretch mode. Either mode can be selectively excited by 
adjusting the shape of the pump pulse through phase 
shaping and/or amplitude shaping of its spectrum. 

Figure ^ depicts the results of a phase-only feedback 
control experiment. Panel a (b) shows the Raman spectra 
produced by the optimal pulse for the symmetric (anti- 
symmetric) stretch mode. This learning search included 
2720 different pulse shapes. Moderate fitness increases 
were observed for either target mode after 25 genera- 
tions. The fitness increase was greater for the symmetric 
mode. This is typical of our searches based on phase-only 
control 18]. 

A Wigner representation of the pulse shape solution 
found by the learning algorithm is plotted in panel c (d). 
The Wigner function is a spectrally resolved field auto- 
correlation: 



W(w,t) = J dui'E(ui - J)E*{u + w')t 



-2iu't 



(4) 



Wigner representations are complete time-frequency 
spectrograms of the optical control field (up to a global 
phase), but the important features leading to control 
of the methanol are obscure. The inability to interpret 
the result easily is typical of many GA search solutions 

dm. 

The principal control analysis of this problem begins 
with a single covariance matrix (Eq. [5J for the entire 
population of pulse shapes evaluated in the two feedback 
experiments. The physical system under control was the 
same in the two problems; only the target was differ- 
ent. We therefore expect the independent searches to be 
nearly spanned by a small number of eigenvector controls. 
The covariance matrix is not simple to interpret because 
the principal components in this problem are widely dis- 
tributed among all of the control settings in the search 
space. However, the essential features of the control 
problem begin to emerge if the covariance matrix is di- 
agonalizcd. This can be seen in Fig. where the eigen- 
values of the covariance matrix are plotted in descending 
order. The algorithm has made large excursions only 
along the few control directions that have large eigen- 
values. The phase functions associated with the three 
largest eigenvalues are shown in Fig. [2fc>. These control 
directions are also the most strongly correlated with fit- 
ness. Figure 13; (d) plots the correlation of eigenvector 
with fitness for the symmetric (antisymmetric) stretch. 
The three most significant eigenvectors are the same as 
the control directions that correlate most strongly with 
the fitness. Therefore we propose that the dimension of 
the search space can be reduced without inhibiting con- 
trol. 

Each pulse shape can now be re-expressed in the eigen- 
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FIG. 2: Panel a: Eigenvalues of the covariance matrix in 
descending order. Panel b: Phase functions associated with 
the three principal control eigenvectors. For the symmetric 
stretch, the principal controls correspond to the third (blue) 
and second (red) eigenvalues. Antisymmetric stretch cor- 
responds to the second (red) and first (green) eigenvalues. 
Panel c (d): Correlation of fitness f with the control vectors 
for the symmetric (antisymmetric) stretch ordered by eigen- 
value. Specifically, the correlation is (< mf > — < rji >< 
f >)/a Vi af. Panel e (f): The correlation in the original basis 
ordered by frequency: (< Sif > — < 8i >< f >)/crj i (T/. 



vector basis. To arrive at the essential features of the 
best pulse, we calculate the projections rjk of the op- 
timal pulse shape onto the principal control directions 
Uk- When these Uk are expressed in the original basis, 
their components are individual discrete frequencies that 
make up the field. Therefore, VjUji 1 ^) produces a 

pulse that contains traits necessary to achieve the target, 
with minimal extraneous features. Figure ^ (f) shows 
the Wigner plot of this essential pulse for the symmet- 
ric (antisymmetric) stretch mode. The essential pulse for 
the symmetric stretch preserves 66% of the original pulse 
shape vector (Xw=i ( w ) = 66%). For the antisymmet- 
ric stretch, this number is 88%. 

This procedure is independent of the specific nature 
of the physical system or the dynamical Hamiltonian; 
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FIG. 3: Top: Magnitude of the Fourier transform of I(t) 
for the optimal pulse shape found by the learning algorithm 
(dashed) and for the essential pulse shape found through prin- 
cipal control analysis (solid) for the symmetric stretch mode. 
Bottom: Similar for the antisymmetric stretch mode. 

however, the principal control directions can be used to 
construct a simplified interaction Hamiltonian, since the 
laser electric field now only depends on a few parameters: 

H(t) = H(E[u j=1 ... k ;t]). (5) 

This can provide important constraints. For example, a 
recent paper demonstrated a control mechanism for SRS 
in methanol based on periodicity in I(t) 01 • The Fourier 



transform of I(t) then reveals the most important Ra- 
man coupling frequencies in the problem 21]. Figure |3 
(top) shows the magnitude of the FT[I(t)] for the opti- 
mal pulse shape found by the learning algorithm for the 
symmetric stretch mode (dashed) and for the same pulse 
shape projected onto the principal control directions u k 
(solid). Figure |21 (bottom) shows similar plots for the 
antisymmetric mode. 

For both modes, projecting the optimal pulses onto the 
principal control directions enhances the frequency com- 
ponents around 3 THz. The coupling frequency found 
through principal control analysis agrees with the model 
based on the mode separation in methanol. Addition- 
ally, comparing the phases of the Fourier transforms for 
each mode yields a phase difference of 7t/2 in the region 
around 3 THz, consistent with the model described in 
reference [l7| . 

In conclusion, we have shown how covariance analy- 
sis of a genetic search algorithm can uncover essential 
features of the dynamical Hamiltonian. Although our 
example involved only phase-shaping of the optical field, 
this technique should be applicable to any system where 
fitness-directed learning algorithms have been used to re- 
veal the path from an initial quantum state to a tar- 
get. Principal control analysis can also be incorporated 
into the experimental search protocol. By discovering 
the principal controls, it should be possible to search the 
space more efficiently, and to test ideas about the sys- 
tem dynamics as the search is proceeding. The method 
should be most useful in cases where the dynamics can 
be described by only a few principal degrees of freedom, 
which are linear combinations of the control parameters 
of the search space. 
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